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ABSTRACT 


This research presents the first accurate three and six Degree of Freedom (DOF) models 
of the small diameter REMUS 100 with cross-tunnel thrusters (CTT). These are the first 
known hydrodynamic models to explicitly consider the impact on navigation of forward 
and aft CTTs. The models presented in this research provide an alternative method to test 
new docking strategies and provide a better understanding of a torpedo-shaped Autonomous 
Underwater Vehicle’s (AUV) dynamics. 


Future AUV missions will have AUVs operating in cluttered and dynamic areas using 
CTTs. This sort of operating environment requires torpedo-shaped vehicles, such as the 
REMUS 100, to maneuver efficiently at slow speeds. The variable-speed CTT models 
in this thesis were originally developed to improve the understanding of vehicle control 
during docking missions. Sub-sea docking stations provide AUVs the ability to remain 
on station for extended periods of time. The extended missions allowed by an AUV 
with docking capabilities greatly improve and increase AUVs’ civilian and military 


applications by decreasing operational costs and increasing on-station time. 


This thesis provides the first variable speed models with CTTs that can be used to develop 
and critique new docking strategies. Installing and extracting a docking system for testing 
is time consuming and expensive. The models created are valuable because they save time 
and resources. These models also provide a cost-effective method for validating control 
strategies and vehicle dynamics. This thesis provides experimental verification of the new 


models and a discussion of the models’ capabilities and limitations. 
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CHAPTER 1: 
Introduction 





1.1 Motivation 

Autonomous Underwater Vehicles (AUV) represent a paradigm shift in U.S. Naval under- 
water operations. The world’s oceans represent an extreme environment that makes the 
use of autonomous systems highly enticing for all actors that strive to perform operations 
in its depths. Rising autonomous systems such as the Remote Environmental Measuring 
Unit (REMUS) 100 vehicle represents not only a feat of engineering prowess but a modern 
platform that can further human exploration and warfare. By furthering the capabilities of 
AUVs, such as the REMUS 100, the Naval Postgraduate School (NPS) furthers the scope 


of modern science and engineering, as well as the Navy’s maritime capacity. 


Work on the REMUS 100 vehicle and its predecessors has been ongoing at civilian and 
military institutions around the globe for the past 20 years. There has been a multitude of 
research into the capabilities of the REMUS vehicle and in producing a simulated model 
that mimics the vehicle’s motion in the underwater environment with high fidelity. Docking 
a REMUS vehicle safely and effectively has been a coveted goal for all parties who operate 
these vehicles. Docking an AUV provides longer on station times and reduces human 
interaction with the vehicle, which in turn increases mission effectiveness and operational 
costs. Research conducted on previous REMUS variants, which had no cross tunnel 
thruster (CTT) capabilities, proved that the vehicle lacked the controllability at low speeds 
to safely and effectively dock the vehicle in a subsea environment. Kongsberg, REMUS’s 
production company, has released a REMUS variant that now implements CTTs in the 
lateral and vertical directions to sustain control of the REMUS vehicle at low speeds. A 
picture of the NPS Center for Autonomous Vehicle Research (CAVR) REMUS 100 with 
CTTs can be seen in Figure 1.1. 





Figure 1.1. Naval Postgraduate School Center for Autonomous Vehicle Re- 
search (CAVR) REMUS 100 Vehicle with CTTs. Source: [1]. 


The REMUS vehicles used in the CAVR have two vertical and two horizontal CTTs. The 
four CTTs generate force on the REMUS using a small propeller at one end of their 
associated tunnels shown in Figure 1.2. The propellers used by the CTTs can operate in 
both directions, however there is a dominant direction of operation where the generated 
force is slightly greater than in the opposite direction. 


The addition of vertical and lateral CTTs to the physical vehicle presents an opportunity for 
new hydrodynamic modeling for small diameter AUVs. Previous models did not provide 
a high-fidelity model of the vehicle’s motion during variable speed operations. The reason 
for the lack of a variable speed model of the REMUS vehicle was twofold: 


1. In previous REMUS variants (no CTTs), the vehicle needed to maintain higher speeds 
to maintain control of the vehicle’s pose. At higher speeds, the control surfaces of 
the REMUS vehicle are effective; allowing for vehicle controllability and system 
modeling. 

2. The software available for the development of a high-fidelity REMUS model that 
handled variable speeds, especially lower speeds where the vehicle is difficult to 
control and hydrodynamic behavior is non-linear, was not available during the creation 
of the original simulated REMUS MATLAB models. 





re] —— , ree 


Figure 1.2. Location of CTTs on a REMUS Vehicle. Source: [1]. 
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This research presents new models of the REMUS vehicle that simulate it’s behavior at lower 
speeds while utilizing CTTs. These models provide the ability to simulate a REMUS vehicle 
performing low speed surveys, CTT maneuvering in confined areas, and hover missions at 
optimal depths for tactical/oceanographic sonar data collection. Also, the models presented 
in this paper bring the goal of reliable autonomous docking of a REMUS AUV closer to 
realization. Using these models to simulate missions prior to deploying the REMUS vehicle 


will provide insight to the vehicle’s capabilities, limitations, and efficient utilization. 


1.2 Problem Description 

Undersea docking stations provide AUVs the ability to remain on station for extended periods 
of time. Extended AUV missions greatly improve and increase widespread applications for 
these vehicles. However, docking stations are not currently in wide spread use due to the 
challenges of consistently and safely docking an AUV. Before consistent docking can occur, 
an accurate hydrodynamic model that can simulate an AUV operating at various speeds 
with CTTs should be produced. This model can provide an accurate simulation of a docking 
maneuver and increase the spectrum of control options available to the vehicle. CTTs allow 
an AUV to hover and perform precise movements at slow speeds. Variable speed modeling 
and CTT functionality is investigated and analyzed in this research. This thesis furthers the 
understanding of AUV hydrodynamic behavior during future docking missions. The added 
capability provided by CTTs will insure damage to both the AUV and the docking station 


can be avoided. 


Several docking missions were attempted in Monterey Bay by the Center for Autonomous 
Vehicle Research (CAVR). These missions used a REMUS variant that did not have a CTT 
capability enabled. The REMUS vehicle was required to maintain a 2-3 knot speed when 
docking to ensure control over the vehicle’s dynamics. At the required speed, the docking 
attempts caused damage to both the vehicle and the docking station, seen in Figure 1.3, and 
were largely unsuccessful. With the results of these missions, it was determined that a CTT 
capability and an accurate hydrodynamic model of the REMUS vehicle with this new CTT 
capability could be used to more reliably dock a REMUS vehicle. The new models provide 
a more concrete understanding of the capabilities of a REMUS vehicle with CTTs. They 
can also be used to better examine the impact of the environment to the vehicle’s behavior 


and help in determining the best control strategy for a docking mission. 
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Figure 1.3. Damage Caused by Unsuccessful Docking Attempts with the 
REMUS at 2-3 Knots. Damage Indicated by Red Circles. Source: [1]. 


To more reliably dock an AUV with CTTs, accurate hydrodynamic models were developed. 
These models allow for accurate simulation and control of the REMUS vehicle. Accurate 
models can be used to help determine optimal trajectories and vehicle limitations under 


various oceanographic conditions. The models presented in this research are: 


1. a three degree of freedom (3DOF) model which represents a REMUS vehicle’s 
behavior when operating two vertical CTTs while at various speeds, 
2. a six degree of freedom (6DOF) model which simulates a vehicle with vertical and 


horizontal CTTs while at various speeds. 


This thesis also provides experimental verification of the new 3DOF and 6DOF models and 
a discussion of their capabilities and limitations. It is distinctive from previous research be- 
cause, while previous work on the REMUS presented hydrodynamic models of the vehicle, 
these models could not accurately reflect the REMUS vehicle’s hydrodynamic behaviors 
while operating at variable velocities and while utilizing CTTs. The new hydrodynamic 
model presented in this thesis is important for accurately representing a vehicle transitioning 


into a hover mode as it slowly approaches a docking station. 


The current docking techniques have an AUV approaching the terminal objective using an 


Ultra Short Baseline (USBL) transponder to determine range and bearing to an acoustic 
transponder mounted on the docking station. Due to a lack of attitude data and low 
periodicity (approximately 1Hz), the USBL does not provide the required resolution of 
position at close range which are necessary for the final portion of a docking maneuver. 
Docking performance could be improved if control algorithms utilize different sensors and 
actuators for precise alignment and positioning at close range, namely a forward-looking 
camera and CTTs. At close range, the REMUS’s control sensor could transition from USBL 
to a forward-looking camera for final alignment and positioning. A combination of the main 
aft thruster and CTTs may be used to position the AUV in front of the docking station. An 
accurate hydrodynamic model that can represent the vehicle’s motion at various speeds and 
with CTTs is needed to develop control algorithms that allow a REMUS vehicle to execute 


a precise, slow-speed approach to a docking station. 


The 3DOF model presented in this thesis focuses on the vertical plane. This research inves- 
tigates the relationships between hydrodynamic coefficients and control surface forces with 
the vehicle’s forward velocity, and the relationship between propeller forces and propeller 
speed (RPM). By coalescing previous research on AUV behavior into 3DOF and 6DOF 
variable speed models, high-fidelity simulations of a REMUS AUV were demonstrated. 


1.3 Literature Review 

Work on AUV hydrodynamic models has been ongoing at civilian and military institutions 
around the globe for over 20 years [2]. There has been a multitude of research into produc- 
ing a model that simulates the REMUS vehicle’s motion in the underwater environment. 


Previous research developed and investigated the following topics: 


. hydrodynamic models that handle single speed dynamics of a REMUS AUV, 
. reduction of cross tunnel thruster efficiency as vehicle speed increases, 
. identification and verification of an AUVs’ hydrodynamic coefficients, 


. accurate simulation of the forces generated by an AUV’s propeller at various RPMs, 


aA BB WO NO eK 


. accurate simulation of the forces on the control surfaces of an AUV at certain speeds. 


1.3.1 Hydrodynamic Models 

There have been several researchers that have developed models of undersea vehicles. [2] 
is a known expert in this field of study and his work developed an accurate and widely 
accepted six degree of freedom model for undersea vehicles. Many of the equations used to 
develop the three degree of freedom model presented in this paper are simplified equations 
from the six degree of freedom model developed by [2]. CTT research has been conducted 
by [1], [3], and [4]. AUV modeling research was pursued by [4], [5], and [6]. All these 
researchers and their work have contributed or influenced the models presented in this 


thesis. 


The hydrodynamic coefficients used in this thesis are adapted from the coefficients developed 
by [7] and the coefficients that were further developed by other researchers. [8] created 
models for a REMUS 100 vehicle in a standard configuration, a REMUS with no CTTs. 
In [1] and [9], hydrodynamic coefficients were further tuned to more accurately reflect 
a REMUS configuration used by the CAVR. The CAVR owns REMUS vehicles that are 
configured to utilize cross tunnel thrusters and house a USBL in the end cap of the vehicle. 
AUV thruster behavior and controllability for a three-state thruster model was analyzed 
by [10] and applied to the REMUS 100 by [8], and then verified by [1]. Further work on 
reduced cross tunnel thruster efficiency as an AUV speed increased was published by [3], 
which was built upon an earlier study by [11], who developed a cross tunnel thruster 


behavioral model at zero speed of advance in 2007. 


None of the prior AUV hydrodynamic models accurately simulated REMUS behavior at 
various speeds while utilizing CTTs. The aim of this thesis was to build off the previous 
work conducted by these researchers and create a hydrodynamic model that could accurately 
reflect an AUV’s behavior in the vertical channel at various speeds. The 3DOF and 6DOF 
models presented account for reduced control surface efficiency at low speeds and build 
an accurate representation of a REMUS AUV’s behavior while operating at slow/variable 
speeds with CTTs. 


1.4 Thesis Organization 
The thesis is organized by the following chapters: 


1. Developing the 3DOF Model for a REMUS - Provides the mathematical derivation 
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of the equations of motion for a 3DOF Model. 

. Programming the REMUS 3DOF Model- An overview of the programming required 
to implement the 3DOF equations of motion for constructing an accurate model of a 
REMUS vehicle. 

. Verification of the 3DOF Model - Results from experiments that were conducted in 
Monterey Bay, CA used to test and validate the 3DOF Model. 

. Developing the 6DOF Model for a REMUS- Provides the mathematical derivation of 
the equations of motion for a 6DOF Model. 

. Programming the REMUS 6DOF Model - An overview of the programming required 
to implement the 6DOF equations of motion for constructing an accurate and more 
detailed model of a REMUS vehicle. This chapter shows the added complexity of 
modeling extra degrees of freedom. 

. Verification of the (DOF Model - Results from experiments that were conducted in 
Monterey Bay, CA used to test and validate the DOF Model. 

. REMUS Controller Design - A discussion of the controllers used in the 3DOF and 
6DOF models and the considerations that went into their design. 

. Behavior and Significance of the Coefficient of Drag (Cd) in a Variable Speed Hydro- 
dynamic Model - Provides insight and modeling considerations specifically regarding 
the impact of slow speed operations 

. Conclusion and Future Work- An exploration of the future work and aspirations that 
the 3DOF and 6DOF models have revealed. 
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CHAPTER 2 
3DOF Model Development 





The focus of this chapter is to develop a vertical channel 3DOF model of a REMUS 100 
AUV with CTTs enabled. A vertical channel model can simulate the vehicle’s control in 
the Z-direction. Due to the constant buoyancy force acting on the vehicle, CTTs must be 
utilized to maintain depth control at slow speeds. The model developed provides insight 
into the behavior of vehicle at slow speeds and illuminates possible challenges to depth 


control while maneuvering a vehicle into a docking station. 


A six degree of freedom (6DOF) model of the REMUS 100 was originally developed by [7] 
and refined by [8]. The 3DOF model in this paper is a simplified version of these 6DOF 
models. However, the previous 6DOF models did not simulate a REMUS 100 vehicle at 
various speeds and utilizing cross-tunnel thrusters. The 3DOF model developed in this 
thesis simulates the behavior of the REMUS 100 vehicle at various speeds and utilizing 
cross-tunnel thrusters. All REMUS 100 mathematical models have been developed using 


AUV equations of motion and kinematics developed by [2] and [12]. 


2.1 3DOF Model Design Considerations 


The aim of this paper is to create a "useful" variable speed model of the new REMUS 100 
variant that utilizes cross-tunnel thrusters. Even when limiting the model to heave, surge, 
and pitch, the 3DOF model is able to show critical vehicle behavior such as the speed at 
which cross-tunnel thrusters need to be utilized to maintain vehicle depth. The 3DOF model 
significantly reduces processing time to run a simulated mission allowing for faster tuning 
and comparison with data collected in the field. Also, a 3DOF model can accurately reflect 
situations in which the REMUS vehicle maintains its depth/altitude when implementing 
cross-tunnel thrusters. The model therefore allows for fast and accurate assessment of 
possible docking missions. External forces in the Z-direction are usually the most limiting 
factor for a successful autonomous docking of an AUV. The model provides the first step 


toward development of a variable speed cross-tunnel thruster 6DOF model. 


The following assumtions for the 3DOF model were adapted from [7] and [8]: 


. No external current is acting on the vehicle. The modeled vehicle is operating in a 
benign operating environment. This assures that the only forces acting on the vehicle 
are from itself and the fluid interaction with the water. 

. No vehicle interaction with waves, external structures, or the sea floor is occurring. 
This assumption further defines benign operating environment and acknowledges that 
any experiments conducted in the field may have these factors acting on the vehicle. 

. Hydrodynamic coefficients of the model change only due to changes in the vehi- 
cle’s velocity. Hydrodynamic coefficients may change due to salinity, pressure, or 
temperature of the water. This assumption is again enforcing the benign operating 
environment for modeling. 

. The vehicle’s center of gravity is aligned with the center of buoyancy when the vehicle 
is at trim. This may not be the case with the real world REMUS, especially if the 
payload and ballasting of the vehicle changes. Often times the center of gravity is 
offset from the center of buoyancy during normal operations. This offset will induce a 
moment in the vehicle that will effect its pitch behavior. Also, if the center of gravity 
is offset from the center line, the vehicles roll behavior may be effected. 

. The vehicle has a slight positive buoyancy. Prior to any operation, the AUV is checked 
to ensure a slight positive buoyancy. The inherent positive buoyancy ensures that if 
the vehicle should loss power while submerged, it will return to the surface so that it 
may be recovered. 

. The vehicle is a rigid body with constant mass. This assumes that the flex inherent 
to the vehicle’s construction and length has a minimal impact on the behavior and 
control of the vehicle. 

. The rudder and dive fins will not stall at any angle of attack. If the rudder and dive fins 
stall, the force equations used to describe forces acting on these control surfaces are 
no longer valid. A stall should not occur for the normal spectrum of vehicle operating 
velocities and the maximum angle of attack for these control surfaces. 

. The assuption that speed of the vehicle does not impact CTT efficiency is necessary 
due to the lack of information on how speed impacts a REMUS’s CTTs. However, 
there has been research that has proven that vehicle velocity does impact CTT effi- 
ciency [3]. The force equations used in this research’s models were developed by [1]. 
The CTT equations for the REMUS were created by measuring the force output of 


CTTs at various RPMs. These force measurements were recorded on a near stationary 
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vehicle. To use these equations to provide a direct force input to the equations of 
motion, the vehicle’s velocity must be assumed to have a negligible impact on CTT 
efficiency. 

9. The forces acting on the vehicle are limited to inertial, gravitational, hydrostatic, 
propulsive, lift, CTT, and hydrodynamic forces. This is a catch-all assumption for 
eliminating all other extraneous forces from impacting the model. 

10. Coupled components of the equations of motion from sway, roll, and yaw are negli- 
gible in the 3DOF model. This allows for the simplifying of the REMUS vehicle’s 
equations of motion and converts a 6DOF REMUS vehicle into a 3DOF REMUS 
vehicle simulation. In the full 6(DOF model presented in Chapter 5, this assumption 


is no longer valid. 


These assumptions provide a standard test setting where environmental and external factors, 


other than hydrodynamic, will have a negligible impact on the model. 
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2.2 Reference Frame and Degrees of Freedom 


r Body-fixed frame 
G 





3 


Figure 2.1. Reference Frame and Degrees of Freedom of a REMUS AUV. 
Source: [8]. 


Figure 2.1 shows the reference frames normally utilized when modeling an AUV. When 
representing the position, orientation, and motion of an AUV, two orthogonal reference 
frames are required. The first reference frame is the global-inertial frame which is a fixed 
frame with respect to the center of the earth. The second reference frame is a body-fixed 
frame with an origin fixed at the vehicle’s center of buoyancy with the principle axes aligned 
with vehicle motion. This body reference frame allows for a more intuitive derivation of 
equations of motion. Both reference frames are right hand reference frames oriented in the 


North, East, and Down directions. 
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2.3 Equations of Motion 
The equations of motion shown in this paper are provided for reference and to show the 
different components of each force and moment considered in the REMUS 100 3DOF 


model. For an in-depth derivation of these equations, please reference [2] and [12]. 


The equations derived in this section are developed in the body reference frame using 
Newton’s second law and Euler’s rigid body equations: 


» F=m*a (2.1) 


SM ale (2.2) 


These forces and moments are a summation of hydrostatic (HS), added mass (A), hydrody- 
namic (HD), lift (L), CTT (CTT) and propeller (P) components. 


»)F = Fus + Fa + Fup + Fy + Ferr + Fp (2.3) 


>) M = Mus + Ma + Muy + Mi + Mcrr + Mp (2.4) 


With these equations, forces and moments for a 3DOF and 6DOF model of a AUV with 
CTTs can be derived. 


Surge Equation 


Table 2.1. Surge Equation Hydrodynamic Coefficients 

















Coefficients | Description Value Units 
Xulul Axial Flow Drag Coefficient | —12.4759 Kg 
Xi Added Mass Coefficient —0.930 kg 
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The following equations derive the summation of forces acting in the X-direction on the 
vehicle. These forces are a summation of hydrostatic, added mass, hydrodynamic, and 


propeller forces. 


> ee Ce Oe Eee (2.5) 


> X =—(W — B)sin(9) + Xjujuluel + Xqt + (1 — t))Tajnjnln| (2.6) 


The propeller force is simplified to a single variable, Fp;opetier, which is a function of 
propeller RPM. 


a a Tp )Tynrln| = F propeller (2.7) 


The alternate force equation is now: 
DX = —(W = B)sin(9) + Xujujeelul + Xatt + Feropetter (2.8) 


Heave Equation 


Table 2.2. Heave Equation Hydrodynamic Coefficients 




















Coefficients | Description Value Units 
Livy Body Lift Force and Fin Lift Coefficient | —28.6 tg 
Zw\w| Cross Flow Drag Coefficient —185.621 Kg 
Zea Fin Lift Force Coefficient —21.37 Kg 5 




















The following equations derive the summation of forces acting in the Z-direction on the 
vehicle. These forces are a summation of hydrostatic, added mass, hydrodynamic, lift, and 
CTT forces. 
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Z= ZHS + ZA + ZHD + ZL + ZCTT (2.9) 


> Z =(W — B)cos 6+ ZyW + Zy|y,W|w| + Zywuw 


+ Zuus,USs + Fyert-CTT, fwd + Fyert-CTT,aft (2.10) 


It is assumed that the two vertical CTTs can be treated as a single force acting in the 
Z-Direction: 


Fyert-CTT, fwd + Fyert-CTT,aft = Fyert-CTT.Total (2.1 1) 


The modified heave equation is now: 


DZ = (W = B) cos 0 + Ziti + Zui wlw] + ZuUw + Zuus, 0s + Frerr—crrroral (2.12) 


Pitch Equation 


Table 2.3. Pitch Equation Hydrodynamic Coefficients 



































Coefficients | Description Value Units 
My Added Mass Coefficient —4.16 kg-m 
Mw Body and Fin Lift and Munk Moment Coefficient | 24.0 kg 
Mug Add Mass Cross Term and Fin Lift Coefficient —10.00 kgm 
My Cross Flow Drag Coefficient 4.00357 kg 
Mus, Fin Lift Moment Coefficient —22;3855 fg 
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The following equations derive the summation of moments contributing to the pitch of the 
vehicle. These moments are a summation of hydrostatic, added mass, hydrodynamic, lift, 


and cross-tunnel thruster moments. 


>) M = Mus + Ma + Map + Mz + Mert (2.13) 


M = —zGW sind + MywW + M,yuw + Mygug + My\y\w|w| 


a Mus, uds as Tyert-CTT,aft + Tyert-CTT, fwd (2.14) 


The vertical thrusters are assumed to be equal distance from the center of buoyancy and are 


assumed to be receiving the same RPM command. Therefore, 


Tyert-CTT,aft = ~Lyert-CTT, fwd (2.15) 


This assumption was made to create a more simplified model of the REMUS 100 which 
still provides an accurate representation of the vehicle’s motion in the Z-direction and pitch 
behavior. The pitch controller in the actual REMUS vehicle ensures that the moments acting 
on the REMUS from the VCTTs are near equal in magnitude. Assuming Tyer;—crr,afr and 
Tyert-CTT,fwd are equal in magnitude and opposite in direction while being equal distance 
from the center of buoyancy essentially provides the same behavior that is created by the pitch 
controller on the vehicle. This assumption reduces the programming and computational load 
of implementing a pitch controller in the model while providing an accurate representation 


of the vehicle’s behavior. 
The total vertical thruster torque is equal to a linear combination of the forward and aft 
vertical thrusters. 


Lert CTE Total = Tyert-CTT,aft + Tyert-CTT, fwd =0 (2.16) 
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This results in the final moment equation: 


dM = -zgWsind + Myw + M,yuw + Myguq + My|njwlw| + Muus,u25s (2.17) 


2.4 Variables Used in the 3DOF Model 


Table 2.4 provides the vehicle specific and environmental values required for coefficient 









































calculations. 

Table 2.4. Critical Values for Coefficient Calculations 
Variable Value Source 
Weight (W) 52.5 kg Bermudez [1] 
Buoyancy (B) 55.8362 N Bermudez [1] 
Length (L) 2.26 m Bermudez [1] 
Hull Frontal Area (A ;) 2.85 « 10°? m? Prestero [7] 
Hull Projected Area (A,) | 0.430512 m? Taylor 
Sea Water Density(p) 1030 kg/m? White [13] 
Sea Water Viscosity (42)! | 0.00141 kg/(m-s) | White [13] 
Diameter of REMUS (D) | 0.190492 m Prestero [7] 
Fin Taper Ratio (t) 0.654 Whicker-Fellner [14] 
Fin Platform Area (Sin) | 6.65 * 10-3 m? Prestero [7] 
Moment Arm (Xfinpost) ?_ | —1.0475 m Taylor 
Fin Lift Slope (c,q@) 3.12 Prestero [7] 














2.5 Vehicle Hydrodynamics, Added Mass, and Thrust Co- 


efficients 
Through analysis of the 3DOF equations in regard to real world vehicle data, it was deter- 


mined that the following coefficients are the major contributors to the equations of motion 
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in a variable speed cross-tunnel thruster model: 


. Xy\u\-Axial Flow Coefficient, 

~ Thnj-REMUS Non-linear Thruster Coefficient, 
. Zy|w|-Cross Flow Coefficient, 

M,y|\w|-Cross Flow Coefficient, 

. Zyy6,- Fin Lift Coefficient, 

» Mius,-Fin Lift Coefficient. 


The following sections provide the calculations and justifications for the 3DOF model’s 


major coefficients. 


The following equation defines the axial drag coefficient X,,,,,: 


1 
Xulul = ~ 5 PCaAy (2.18) 
p and Af are both constants, however the coefficient of drag (Cz) will change with the 
vehicle’s velocity depending on the Reynolds number. To better interpret how this Cy will 
be affected, the Reynolds number in the X-direction for the REMUS needs to be calculated 


over its entire operating regime [13]: 


_ p|U|D 
L 


Re 





(2.19) 


Diameter (D) was used instead of the vehicle cross-sectional area to ensure that the smallest 
Reynolds number that could be seen by the vehicle was calculated. This calculation 
provided the most conservative analysis of the fluid dynamic properties of the REMUS 


vehicle operating at variable speeds. 


The value of the Reynolds number determines turbulent or laminar flow by the relationships 
seen in Table 2.5. 


As can be seen in Figures 2.2 and 2.3, nearly the entire range of operating velocities for 
the REMUS vehicle is in the turbulent region. The REMUS vehicle transitions from the 
turbulent region at approximately 0.071863 m/s. Therefore, the 3DOF model assumes there 
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is slight to moderate Reynolds number dependence when calculating Cg for the entire range 


of velocities. 


4 
ae 


a — 
(o>) © (=) NO 
1 1 


Reynolds Number (Re) 
_ 











0 l l l 
0 0.2 0.4 0.6 0.8 1 
Velocity (m/s) 


Figure 2.2. Reynolds Number vs. REMUS velocity 


Table 2.5. Reynolds Number Regimes. Adapted from [2], [13]. 





Re Range Behavior 

0<Re<l Highly viscous laminar creeping motion 

1 < Re < 100 Laminar, strong Reynolds number dependence 
100 < Re < 10° | Laminar, boundary layer theory useful 

10° < Re < 10° | Transition to turbulence 

10* < Re < 10° | Turbulent, moderate Reynolds number dependence 
10° < Re < co | Turbulent, slight Reynolds number dependence 
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Figure 2.3. Transition Thresholds of the Reynolds Number for the REMUS 


Due to the Reynolds number behavior for the REMUS vehicle, seen in Figures 2.2 and 2.3, 
Cg will remain constant at all speeds for the 3DOF model. Cy for this research was taken 
from previous research conducted by [13] and [15]. For a blunt spherical object such as the 
nose cone of a REMUS, with a Re > 10* and with D/L ~ 10 (REMUS D/L = 11.59), 
Cq ~ 0.82 according to [13] or Cg ~ 0.85 according to [15]. Cg = 0.85 was used in the 


model presented in this paper. 


Using the Cg above, X,,,, can be calculated for the entire operating region of a REMUS 


vehicle which implements cross-tunnel thrusters. 


1 
Xuiu) = — #103082. « 85 + 2.85 « 102m? (2.20) 
2 m 
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k 
Xu) = -12.4759—2 (2.21) 
m 


The cross flow drag coefficient (Z,,|,,|) can be broken up into two components, Zy|w|,Body 
and Zy\w|,Fins,s eXpressed in equation 2.22: 


LZw\w| ™ Zw\w|,Body + Zw\w|,Fins (2.22) 


The Reynolds Numbers in the Z-direction will be much higher than the Reynolds Numbers 
in the X-direction due to the larger associated area of fluid interaction. It was proven in 
the X,,|,; section that the REMUS vehicle is always operating in a turbulent flow regime. 
Since the vehicle is always in a turbulent flow, the cross flow drag coefficient Cg. can be 
determined and the assumption that Re > 10* in the Z-direction is valid. 


For a cylinder with a Re > 10* and with D/L ~ 10 (REMUS D/L = 11.59), the Cac 
is approximately 0.82 according to [13] or Cac is approximately 0.85 according to [15]. 
Cac = 0.82 was used in the model presented in this paper. Giving us the following value 


for Zw|w|,Body: 


1 k 
ZoujBody = ~5PCacAp = -181.8705-4 (2.23) 
Using the cross flow drag coefficient (Cyr) developed by Whicker and Fehlner [14] and the 


fin taper ratio (t) the following equations are derived: 


Cap = 0.1 + 0.7t = 0.1 + 0.7 * 0.654 = 0.56 (2.24) 
1 k 
Zojwhins = 2* 5PCapSfin = ~3.15-—* (2.25) 
m 


Therefore, 
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k 
Zulu) = —185.621— (2.26) 
m 


The cross flow coefficient (Z,,|,,;) above is valid for the entire range of REMUS operating 
velocities since the vehicle is always operating in the turbulent flow regime where there is 


slight to moderate Reynolds number dependence. 


Cross flow drag (M,,|,|) is calculated in a similar way as the cross flow coefficient (Z,,|,): 
Mow} es My\w|,Body + My\w\,Fins (2.27) 
where, My|w|,Body = 9, due to the assumption that the center of buoyancy lies on the center 


of the vehicle, equidistant from the nose and tail. With this assumption, the following 


equation is valid: 


Moy|\w| = Molw|,Fins (2.28) 
where 
Mo|w\,Fins = 2 * Xfinpost * (-5pCas fin) (2.29) 
therefore, 
Myj\w| = 4.00357kg (2.30) 


The cross flow drag coefficient (M,,|,,|) above is valid for the entire range of REMUS 100 
operating velocities since the vehicle is always operating in the turbulent flow regime where 


there is slight Reynolds number dependence. 


The fin lift force (Z,,5;) is determined by the following equation: 


Ze 


1 
Zuués = 2 * (5 pCLoS fin) = —21.37kg/(m- rad) (2.31) 


The fin lift moment (M,,,55) is determined by the following equation: 


1 
Muss =2* (5 PCLaS fin finpost) = —22.3855kg/rad (2.32) 


2.6 Other Coefficients 


The following coefficients are the minor contributors to the equations of motion in a variable 
speed cross-tunnel thruster model. These coefficients were deemed minor due to the lower 


velocities and accelerations associated with the listed coefficients: 


. Xudotr-Ahdded Mass Coefficient, 

. Zwdor-hdded Mass Coefficient, 

. Zyw-Body Lift and Fin Lift Coefficient, 

. M,w-Body Lift and Fin Lift and Munk Moment Coefficient, 
. Myg-Added Mass Cross Term and Fin Lift Coefficient. 


nN BW NY 


The behavior of these minor contributors regarding variable velocity is not addressed in this 
research. However, the constants for these coefficients were derived in previous research 
for the REMUS vehicle and can be found in [1] and [7]. 


2.7 Direct Force and Torque Inputs 
The following are the direct force and torque inputs applied to the equations of motion for 
the 3DOF model: 


1. F Propeller =Xp, 
2.. Fyert=CTT Total = ZCTTs 
3. Tyert-cTT,Total = Mcrr. 


F propeller (Xp) for the REMUS 100 was derived through experimentation utilizing a FUTEK 
strain gauge. The following equation relates vehicle RPM (n) to Fp;opeijer in Newtons: 
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Fpropeller = 2.93 X 107 x n* — 2.69 x 10°? xn? 


+ 7.23 x 107° x n* +0.0105 xn (2.33) 


The equation for Fp; opeiier has the behavior seen in Figure 2.4. 


Fae 2.93e-12*n* + -2 69e-09*n? + 7.23e-06*n2 + 0.0105*n 


> 
o 


Force (N) 
3 
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Figure 2.4. Fpropelier VS. Propeller RPM. Source: [1]. 


Fyert-cTTTotal (Zcrr) for the REMUS 100 was derived through experimentation utilizing a 
FUTEK strain gauge. The strain gauge experiments were conducted to find the force gen- 
erated by one cross-tunnel thruster at a certain RPM (n) in the downward direction (Faown). 
The experiments yielded a different force equation when the thrusters were operating at 


a certain RPM (-n) in the upward direction (F,,). With this research assumptions, the 
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following relation is true: 


2% Fuown, Whenn > 0 
Fyert—-CTT Total = 
—2* Fyp, Wwhenn <0 


The following equations relate vehicle RPM (n) to Faown and Fy, in Newtons: 


Frown = 3.89 X 107 x n* +3.5 x 107!! x n3 


+1.81x 1077 xn? +3.3x107*xn (2.34) 


Fup = 6.85 X 10° xn* +7.3 x 107"! xn? 


+ 9.26 x 10°78 xn? +4.97x 104+ xn (2.35) 


The equations F,,, and Faown have the behavior seen in Figures 2.5 and 2.6. 
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Force (N) 





= -6.85e-15*n4 + 7.3e-11*n? + 9.26e-08"n2 + 0.000497*n 
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Figure 2.5. Fy» vs. Cross-Tunnel Thruster RPM. Source: [1]. 
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Feat 3.89e-15*n* + 3.5e-11*n? + 1.81e-07*n2 + 0.00033*n 


Force (N) 


QO Test Data 
Polyfit 
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RPMs 


Figure 2.6. Faown vs. Cross-Tunnel Thruster RPM. Source: [1]. 


Torque (Tyert—CTT,Total Ot Mcrr) due to output of the horizontal cross-tunnel thrusters in the 
model is assumed to have the following property: 


[Tvert-CTT,aft\ = [Tvert-CTT, fwdl (2.36) 


Therefore, 


Tyert-CTT,Total =0 (2.37) 
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CHAPTER 3: 
Design of a 3DOF Computer Model 





This chapter is an overview of the design of the computer model created in MATLAB 
Simulink for the REMUS 100. The 3DOF model presented in this chapter is partially 
verified through experiments in Monterey Bay, CA. This chapter is included in this thesis to 
provide a more in-depth understanding of the model for future researchers that continue the 
verification and validation process. It will also provide an overview of how an individual 


could create their own variable speed model of an AUV. 


3.1 Equations of Motion 
The first step in creating a variable speed model was implementing the equations of motion, 
found in Chapter 2, into a Simulink 3DOF model. The implementation for the surge, heave, 


and pitch equations are depicted in Figures 3.1, 3.2, and 3.3. 
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Figure 3.1. 3DOF Surge Equation 
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Figure 3.2. 3DOF Heave Equation 


The Bermudez force equation block in Figure 3.2 contains the CTT RPM to Force equation 
derived by [1]. The input to this block is the RPM for both VCTTs and the output is the 
Force (N) that is generated by the simulated VCTTs. 
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Figure 3.3. 3DOF Pitch Equation 


3.2 Transformation From Body Fixed to Global Reference 
Frame 


The next step was transforming the forces in a body fixed coordinate frame into motion in the 
flat earth global reference frame. This transformation was easily performed using a function 
block from the Aerospace Toolbox in Simulink which converted the forces generated by the 
equations of motion into motion in the global frame. The Aerospace function block is seen 
in Figure 3.4. 
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Figure 3.4. Equations of Motion Transformed from Body Fixed to Global 
Reference Frame 


3.3. Implementation of the REMUS Controller 

The third step in the design process required commands (i.e. depth and RPM commands) 
to compare, through the use of a feedback loop, the current position of the model to the 
commanded position. To verify that the model’s behavior was similar to the actual vehicle’s 
behavior, real world commands from previous missions were required. The model was then 
verified by issuing these recorded real world commands to the model and comparing the 
model’s behavior to the recorded behavior of the REMUS vehicle operating in Monterey 
Bay. A detailed analysis of the model’s performance can be seen in Chapter 4. The real 
world commands were imported into the Simulink model as an array and were issued to the 
model at the rate seen by the actual REMUS vehicle. The commands issued to the model’s 
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controller were depth and propeller RPM commands. The implementation of the real world 


commands is shown in Figure 3.5. 
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Figure 3.5. Controller Importing Real World Commands. 


The real world commands were accepted by the model’s controller, seen in Figure 3.5. 
The controller in the model is actually a group of PID controllers, controlling Propeller 
RPM, Depth, and Cross-Tunnel Thruster RPM of the model. The modeled controller 
mimics the PID controllers implemented by the REMUS vehicle, however it is not an exact 
copy of the REMUS vehicle’s controller due to the actual REMUS controllers’ design and 
parameters being proprietary information. However, even with a slightly different controller 
implemented in the model, the actuation of the control surfaces, thrusters, and propeller 
occur so rapidly that the model’s behavior was not affected significantly (refer to model 


performance in Chapter 4). 
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3.4 Complete Model 


The equations of motion, body-fixed to global reference frame transformation, real-world 
command inputs, and PID controllers comprise a full 3DOF model to simulate the AUV 


vertical plane motion at various speeds, depicted in Figure 3.6. 


REMUS Control System 


XZ Graph 








‘3D OF (Body Axes) 


Figure 3.6. Complete Variable Speed REMUS Model Utilizing Cross- Tunnel 
Thrusters. 


The reason the program seen in Figure 3.6 works at various vehicle velocities is that the 
equations of motion are recomputed at each iterative time step, whose duration is set in the 
model’s parameters. This iterative calculation process is a stark change from previously 
presented models of the REMUS vehicle and allows for a more detailed and representative 
analysis of the vehicle’s behavior. The time step used for this model was 0.002 seconds. 
At this time step interval the modeled vehicle can accept all commands from real world 
missions and compute the vehicle’s behavior in real-time or faster than real-time. The 
model’s fixed time step is limited by the data rate of the vehicle’s sensors and commands. 


The 0.002 second time step was selected because at this time interval the model is able to 
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accept these real world data and commands while calculating the vehicle’s behavior to a 


high accuracy. 
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CHAPTER: 
REMUS 3DOF Model Verification 





Presented in this chapter are the verification experiments of the 3DOF model’s depth 
behavior. Depth control at the early stages of this research was considered to be the 
parameter of greatest concern when modeling a REMUS vehicle. Depth control is a 
difficult element to model due to the constant positive buoyancy force acting on the vehicle 
and the changing efficiency of the dive fins at various speeds. To maintain depth without 
CTTs, the REMUS must maintain a certain speed to ensure adequate flow over the dive fins. 
At slow speeds, CTTs are the only way to maintain depth control. Slow speed maneuvering 


is paramount when modeling an AUV attempting to dock. 


4.1 Variable Speed Verification 

The 3DOF model verification was conducted in two parts. The first part of the verification 
process was conducted by obtaining REMUS commands and data from experiments in the 
field. A thousand seconds of depth and RPM commands from a REMUS vehicle operating 
in Monterey Bay were imported into the 3DOF model. In the first part of the verification, 
cross-tunnel thrusters were not enabled on the vehicle. However, depth and RPM commands 
did vary during the analyzed thousand second window. The 3DOF model’s behavior in the 
Z-direction was then plotted against the actual vehicle behavior and a sensitivity analysis 


was conducted. The results are represented in Figures 4.1, 4.2, and 4.3. 
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Figure 4.1. REMUS Estimated Forward Velocity 
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Figure 4.3. REMUS Real World vs. Modeled Behavior 


In Figure 4.3, at the beginning of the plot, there is a noticeable difference between the 
REMUS vehicle’s actual depth and the model’s depth. This difference is due to the initial- 
ization of the model at a depth and speed slightly different than the actual vehicle. As the 
model’s depth controller reacts to the depth commands imported from the actual mission 
and the buoyancy force, the model’s dive fin angle changes and its behavior stabilizes to the 
depth of the actual vehicle. Also, as the model’s velocity more closely follows the actually 


vehicle’s velocity, it is able to decrease the error in its depth control behavior. 


At approximately 840 seconds, the actual vehicle loses some of its depth control capability. 
This behavior is caused by a reduction of flow over the dive fin surface. The reduction of 
flow across the dive fins may have been caused by a difference in current at this depth and a 


decrease in vehicle velocity. The vehicle’s propeller RPM at this time was 799 RPM which 


40 


is low in its normal operating range without CTTs enabled. 


The mission profile for Figures 4.1, 4.2, and 4.3 was selected because the actual REMUS 
vehicle operated at various depths between 25 and 30 meters while also reducing RPM from 
1640 RPM to 790 RPM and then increasing RPM back to 1640 RPM. This provided the 


model a variable range of depths and velocities against which it could be analyzed. 


Through performing a sensitivity analysis, it was determined that the mean of the error 
between the real world REMUS behavior and the model behavior was -12.8 cm with a 
standard deviation of 34.4 cm. This error calculation was for the entire thousand second 


mission. 


4.2 Cross-Tunnel Thruster Verification 

The second part of the 3DOF model verification process tested the cross-tunnel thruster 
portion of the model. Another 1000 second window was selected from field experimentation 
conducted in Monterey Bay. The REMUS vehicle was operating at various speeds with 
varying cross-tunnel thruster inputs. Again, depth and RPM commands were fed to the 
3DOF model and the model’s behavior in the Z-direction was plotted against the actual 
vehicle’s behavior and a sensitivity analysis was conducted. The results are seen in Figures 
4.4, 4.5, 4.6, and 4.7. 
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Figure 4.4. REMUS Real World Cross-Tunnel Thruster Commands 
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Figure 4.5. REMUS Estimated Forward Velocity 
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Figure 4.6. REMUS Real World RPM Commands 
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In Figure 4.7, at the beginning of the plot, there is a noticeable difference between the 
REMUS vehicle’s actual depth and the model’s depth. Many of the initialization considera- 
tions that applied to the CTT disabled mission also apply to this CTT enabled mission. The 
vehicle velocity and flow over the dive fins do not heavily impact the actual and modeled 
vehicles’ depth behavior in a CTT enabled mission. The actual and modeled vehicles are 
able to achieve the commanded depth faster than CTT disabled vehicles. There is a slight 
overshoot of the commanded depth by the modeled REMUS vehicle due to the implementa- 
tion of a less efficient depth controller. This overshoot may be corrected by a closer tuning 
of the depth controller in the model. The model implemented bang/bang control of the 


VCTTs which is an energy inefficient implementation of a controller but closely represented 


the VCTT behavior of the actual REMUS vehicle. 
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Through performing a sensitivity analysis, it was determined that the mean of the error 
between the real world REMUS behavior and the model behavior was 8.8 cm with a 
standard deviation of 43.4 cm. This error calculation was for the entire thousand second 


mission. 


The verification experiments of the model’s depth control proved that the equations and 
controllers utilized in the 3DOF model provided a high fidelity representation of the depth 
behavior of a REMUS vehicle. This model can be used to simulate the M, X, and Z- 
direction behaviors of a REMUS vehicle. It can be useful in determining a successful 
docking strategy and optimal docking trajectory. This model can also aid in the design of 
future docking stations by providing designers an accurate representation of the vehicle’s 
depth control. Docking station capture cages can now be designed to ensure the maximum 
probability of the vehicle successfully docking based on the selected docking strategies and 


trajectories. 
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CHARTER 5: 
6DOF Model Development 





Once a 3DOF model of an AUV was designed, the research was expanded to create a 6DOF 
model for a REMUS vehicle operating at various speeds with CTTs. The 6DOF model 
presented in this chapter includes the equations and coefficients calculated in Chapter 2, but 


also includes equations and coefficients presented by [7] and [8]. 


5.1 Equations of Motion 

The assumptions found in Chapter 2 were also used for the 6DOF model. One added 
assumption for the 6DOF model is that roll of the REMUS vehicle is negligible. This 
assumption reduces the computational burden of one degree of freedom. This assumption 
is valid because the roll of a REMUS vehicle varies between -3 to 3 degrees in a normal 
mission. This is a slight difference in roll angle and the variation does not appreciably effect 
the vehicle’s overall behavior. The reference frames presented in Chapter 2 are also used in 
the development of the 6DOF model. A detailed derivation of the 6DOF model equations 
is provided in [2] and [12]. 


Surge Equation 


Table 5.1. 6DOF Surge Equation Hydrodynamic Coefficients 


























Coefficients | Description Value Units 
Xulul Axial-Flow Drag Coefficient —12.4759 Kg 
Xi Added Mass Coefficient —0.930 kg 
Xwq Added Mass Cross-Term Coefficient | —77.8 fg 
Xaq Added Mass Cross-Term Coefficient | —1.93 kgm 
Ry Added Mass Cross-Term Coefficient | 35.5 Kg 
Xrr Added Mass Cross-Term Coefficient | —1.93 kgm 
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The following equations derive the summation of forces acting in the X-direction on the 
vehicle. These forces are a summation of hydrostatic, added mass, hydrodynamic, and 


propeller forces. 
DX = Xus + Xa + Xuv + Xp (5.1) 
> X = —(W — B)sin(@) + Xjulu| + Xai + Xvgwg 
+ Koad +X,port+ Xr +(1- Tp)Tnjn|n| (5.2) 


The propeller force is simplified to a single variable, Fp,;opetier, which is a function of 
propeller RPM. 


(1 = Tp)Tn|njr|n| = F Propeller (5.3) 


The alternate force equation is now: 


SX = —(W = B)sin() + Xujujuluel + Xatt + Xvgwg 


+ Mag + Xp, Ur + More a F propeller (5.4) 
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Sway Equation 


Table 5.2. 6DOF Sway Equation Hydrodynamic Coefficients 









































Coefficients | Description Value | Units 
Yin Body Lift Force and Fin Lift Coefficient —28.6 Kg 
Yoo Cross-Flow Drag Coefficient —2850 Kg 
Yyr| Cross-Flow Drag Coefficient 0.632 oe 
Y; Added Mass Coefficient -77.8 | kg 
Y;. Added Mass Coefficient 4.16 kgm 
Vie Added Mass Cross-Term and Fin Lift Coefficient | 5.22 fg 
Yuu, Fin Lift Force Coefficient 9.64 a 7 

















The following equations derive the summation of forces acting in the Y-direction on the 
vehicle. These forces are a summation of hydrostatic, added mass, hydrodynamic, lift, and 


horizontal cross-tunnel thruster forces. 


SY = Yas +¥a+ Yao + Yi + Yorr (5.5) 


», Y =(W- B)cosO@sing + Y,,uv + Y,jv|v| + Y,r|r| + Yo + YF 


2 
+ Y,pur + Yypwp am YngP4 + Yuud,U Or + Fhorizontal-CTT, fwd + Fhorizontal-CTT,aft (5.6) 


Since vehicle roll is assumed to be 0, the sway equation can be simplified to the following: 


Dt = Y,puv + Y,,)v|v| + Y4-r|r| + Yo + Yer + Y,-ur 


2 
+ Yuué,U 6, + Fhorizontal-CTT, fwd ad Fhorizontal-CTT,aft (5.7) 
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Heave Equation 


Table 5.3. 6DOF Heave Equation Hydrodynamic Coefficients 









































Coefficients | Description Value Units 
Li Added Mass Coefficient -77.8 kg 
Li Added Mass Coefficient —4.16 kgm 
Lug Added Mass Cross-Term and Fin Lift Coefficient | —12.22 Kg 
ZiAa| Cross Flow Drag Coefficient -0.632 | 
Zuw Body Lift Force and Fin Lift Coefficient —28.6 Kg 
Zw|w| Cross-Flow Drag Coefficient —185.621 Kg 
Ziiué Fin Lift Force Coefficient —21.37 _o 5 








The following equations derive the summation of forces acting in the Z-direction on the 


vehicle. These forces are a summation of hydrostatic, added mass, hydrodynamic, lift, and 


vertical cross-tunnel thruster forces. 


SZ = Zns + Za + Zap + Zi + Zerr 


oy Z =(W — B)cos @cos ¢ + ZywW + Zyq + Zuquq + Zypvp + Z-prp 


(5.8) 


7 Zw|w|W|w| + Zq\q\ (14 + Zywuw + Ziaie Os + Fyert-CTT, fwd + Fyert-CTT,aft (5.9) 


It is assumed that the two vertical cross-tunnel thrusters can be treated as a single force 


acting in the Z-Direction: 


Fyert-CTT, fwd + Fyert-CTT,aft = Feni=CTr Total 


The modified heave equation is now: 
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(5.10) 


ss Z =(W — B)cos 6 cos ¢ + ZywW + Zgq + Zuquq + Zypvp + Z-prp 


+ Zw|w|w|w| f Zal\a\ 4 + Zywuw + Zuus,USs + Fyert-CTTTotal (O.11) 


Since vehicle roll is assumed to be O in this model, the heave equation can be further 


simplified to the following: 


>» Z =(W — B)cos 6 + Zyw + Zyq + Zugug + Zy|y\wl|w| + Zgigialal 


2 
+ Zyyuw + Lund, Os + Fyert-CTT,Total (5.12) 


Roll Equation 
The following equations derive the summation of moments contributing to the roll of the 
vehicle. These moments are a summation of hydrostatic, added mass, hydrodynamic, and 


propeller moments. 


))K = Kus + Ka + Kup + Kp (5.13) 


For a full 6DOF model, the sum of the moments for roll would be defined by the following 


equation: 


>; K =-zgW cos 0 sin 6 + Kpp + Kpip\p|pl + Qnininlnl (5.14) 


However, for a torpedo shaped vehicle such as the REMUS, roll can be assumed to be 
negligible in most circumstances. Therefore, the sum of roll moments for the model 
presented in this research is represented by the following equation: 


SK -O (5.15) 
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Pitch Equation 
Table 5.4. 6DOF Pitch Equation Hydrodynamic Coefficients 


















































Coefficients | Description Value Units 
M, Added Mass Coefficient —30 a 
Maiq Cross-Flow Drag Coefficient —188 Es 
My Added Mass Coefficient —4.16 kg-m 
Mw Body and Fin Lift and Munk Moment Coefficient | 24.0 kg 
Mug Add Mass Cross Term and Fin Lift Coefficient —10.00 kgm 
My Cross-Flow Drag Coefficient 4.00357 kg 
Mus, Fin Lift Moment Coefficient —22.3855 Aa 





The following equations derive the summation of moments contributing to the pitch of the 
vehicle. These moments are a summation of hydrostatic, added mass, hydrodynamic, lift, 


and cross-tunnel thruster moments. 


>) M = Mus + Ma + Muy + Mi + Mcrr (5.16) 


> M = ~zGW sin + Myw + Mgq + M,,uw + M,pop + Myprp 


oe Muguq + Mo|w\v|w| + Maqialal + Muus,U5s + Tyert-CTT,aft + Tyert-CTT, fwd (5.17) 


The vertical thrusters are assumed to be equidistant from the center of buoyancy and are 
also assumed to be receiving the same RPM command. The reasoning for this assumption 


is found in Chapter 2 in the Pitch Equation subsection. Therefore, 


Toert-CTT,aft = ~Lvert-CTT, fwd (5.18) 


The total vertical thruster torque is equal to a linear combination of the forward and aft 
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vertical thruster torques. 


Tert=CTT Toial = Tyert-CTT,aft + Tyert-CTT, fwd =0 (5.19) 


This results in the following moment equation: 


)) M = -<GW sin 8 + Mytd + Mag + Muyuw + Mypop + Mrprp 
+ Mugq + Moy|w\v|v| + Mogi \4 + Myus,U6s (5.20) 


Considering that roll is negligible for the REMUS vehicle, the pitch equation can be further 
simplified to the following: 


DM = -<GW sin 8 + Mytd + Mag + Myyuw 
+ Mugtq + My\w\w|w| + Maqialal + Mus, 6s (5.21) 


Pitch Behavior Parameter 

It was found during the creation and verification of the 6DOF model that the righting 
force that stabilizes pitch of the REMUS 100 was too weak to counteract the impact of the 
cross-term coefficients in the equations of motion. To correct this error, a new coefficient 
was introduced called the pitch behavior parameter, Mpg. This parameter represents the 
vehicle’s inherent tendency to return its pitch to a trim position. This inherent behavior is 
a result of the center of gravity being positioned below the vehicle’s center of buoyancy. 
Whenever the vehicle’s center of gravity is disturbed to a position other than directly under 
the center of buoyancy, the vehicle witnesses a large correcting moment that returns it to 
a trim position. Previous models didn’t adequately simulate this rapid increase in pitching 


moment that forcefully and abruptly positions the vehicle to a trim position. 


One explanation for the use Mpz is that 6DOF equations for pitching may have an error 


aS 


in the —zGW sin @ term. This term only takes into account the force caused by the weight 
of the vehicle (W) and the righting moment arm. However, other studies on underwater 
vehicles show that the buoyancy force (B) has a substantial impact on the righting moment 
of an underwater vehicle, illustrated in Figures 5.1 and 5.2. 


Buoyant Force 





Weight 


Figure 5.1. Center of Gravity (CG) and Center of Buoyancy (CB) in Equi- 
librium (Vehicle at Trim). Source: [16]. 
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Buoyant Force 


Figure 5.2. Weight and Buoyancy Working Together to Return the Vehicle 
to Trim. Source: [16]. 


If buoyancy force (B) does have the effect that is expressed in Figure 5.2, the —zGW sin 6@ 
term should be changed to —zg(W + B)sin@. This will more than double the righting 
moment caused by pitch of vehicle (6). Also, there is speculation that the distance between 
center of gravity and center of buoyancy (zg) is four to five times larger than the 1.96 cm 
defined by previous models. The largest zg possible is the REMUS vehicle’s diameter, 
19.05 cm. When considering a maximum 2g and an influence of buoyancy force, there is 
the potential of increasing the magnitude of righting moment up to 18 times the original 


value. However, even with maximum Zg and the addition of buoyancy force to the model, 
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the model still does not have the righting moment behavior that is witnessed in the real 
world. There is possible influence from the long length and construction of the REMUS 
vehicle, as well as its distribution of ballast along this entire length, that creates an even 
stronger righting moment. The proposed moment equation with the buoyancy force taking 


part in the righting moment is: 


Se M = -zo(W + B)siné + Myw + Mgq + M,puw 
+ Mygtg + My|y\w|w| + My\q\4141 + Myus,U?5s (5.22) 


A simple experiment only needs to be performed to comprehend and grasp the importance 
of this trim correcting moment. If the REMUS were placed submerged in a pool and shoved 
in any direction or to any angle of pitch, the vehicle will inherently and quickly return to 
its stable trim position. The REMUS was specifically designed so that it would return to 
trim without external forces during a disturbance. Mpg was developed and implemented 
to model this inherent and rapid righting moment. This research intended to utilize the 
equations of motion developed and adhered to by previous researchers, thus a gain of Mpg 
was introduced to the pitch equation. The influence of Mpg will greatly vary depending 
on the payload and the ballast of the vehicle. Based on behavior witnessed in experimental 
tests and through verification with the 6DOF model, the NPS REMUS vehicles have a 
Mpg ~ 50. Mppz = 50 was utilized in the 3DOF and 6DOF models and the parameter can 
be seen in the corrected pitch equation, Equation 5.23. Further study and analysis of the 
pitch equation and pitch behavior is required to fully correct the 6DOF model. 


yim = —MppzGW sind + Myw + Mgq 


+ M,yuw + Mugug + My|w\w|w| + Maigialal + Myus, Us (5.23) 
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Yaw Equation 


Table 5.5. 6DOF Yaw Equation Hydrodynamic Coefficients 


















































Coefficients | Description Value Units 
Ni Added Mass Coefficient 4.16 kg-m 
N;: Added Mass Coefficient 4.88 ee 
Nw Body and Fin Lift and Munk Moment Coefficient | —24 kg 
Nur Added Mass Cross-Term and Fin Lift Coefficient | —2 kgm 
Nolo | Cross-Flow Drag Coefficient —3.18 kg 
Ny\r| Cross-Flow Drag Coefficient —245 am 
Nuwué, Fin Lift Moment Coefficient —22.3855 Ag 








The following equations derive the summation of moments contributing to the yaw of the 
vehicle. These moments are a summation of hydrostatic, added mass, hydrodynamic, lift, 


and cross-tunnel thruster moments. 


SIN = Nus + Na+ Nav + Nz + Nerr (5.24) 


>, N = Nj0 + Nef + Nyuv + Nypwp + Npgpg + Nurur 


2 
+ Nyjo|P2| “As Ny |r Ir | + Nuns, 6, + Thorizontal-CTT,aft + Thorizontal-CTT, fwd (5.25) 


Finally, since roll is negligible for the REMUS vehicle, the yaw equation can be further 
simplified to the following: 


>» N =N,b + Ne? + N,,uv +N, ur 


2 
an Nolo |0| + Nyr\F IF | + Nuns, 6, + Thorizontal-CTT,aft + Thorizontal-CTT, fwd (5.26) 
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Vehicle Hydrodynamics, Added Mass, and Thrust Coefficients 

Chapter 2 provides a derivation of the specific coefficients needed to model the REMUS 
vehicles owned by the CAVR. These are the vehicles that are used for the verification of the 
6DOF model. For a more detailed derivation of the coefficients used in the 6DOF model, 
please refer to [1], [7], and [8]. A list of all coefficients used in both the 3DOF and 6DOF 


models are provided in the appendix. 


5.2 Conclusion 

The 6DOF model is inherently more complex than the 3DOF model due to the cross- 
coupling terms in the equations of motion. This research started with the development of a 
3DOF model to ensure that each term in these equations were understood and correct values 
for the coefficients were used. Once the equations were analyzed and the correct values 
were selected, a 6DOF model that accurately represented the vehicles used by the CAVR 
was created. The six equations represented in this model are highly impacted by each other. 
Changing a cross-coupling coefficient in one equation may have large residual effects in 
others. Each manipulation of an equation or coefficient must be critiqued in a manner that 


takes into account the changes that will occur to the model as a whole. 
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CHAPTER 6: 
6DOF Computer Model 





The 6DOF model presented in this chapter is an extension of the 3DOF model in Chapter 
3. However, there are some key differences and increased sophistication in the 6DOF 
model. The graphics and discussion are included in this thesis to provide a more in-depth 
understanding of the 6DOF model for future researchers and will be of assistance to any 


verification process that is to be conducted on the 6DOF model. 


6.1 Equations of Motion 

As with the 3DOF model, the equations of motion are the backbone of the programming 
process. Understanding the cross-term hydrodynamic coefficients becomes increasingly 
important when scaling the model from a 3DOF to a 6DOF model. The implementation of 
the surge, sway, heave, roll, pitch, and yaw equations are seen in Figures 6.1, 6.2, 6.3, 6.4, 
6.5, and 6.6. 
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Figure 6.1. 6DOF Surge Equation 
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Figure 6.2. 6DOF Sway Equation 
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Figure 6.3. 6DOF Heave Equation 


Figure 6.4. 6DOF Roll Equation 
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Figure 6.5. 6DOF Pitch Equation 
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Figure 6.6. 6DOF Yaw Equation 


6.2 Model Transformations 

Similar to the 3DOF model, the forces acting on the body had to be converted to positions, 
velocities, and accelerations in the flat earth global reference frame. There is a 6DOF 
function block in the MATLAB Simulink Aerospace Toolbox that transforms the outputs 
of the equations of motion to positions, velocities, and accelerations in the flat earth global 
reference frame. The inputs and outputs associated with the 6DOF function block are seen 


in Figure 6.7. 
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Figure 6.7. 6DOF Equations of Motion Transformed from Body Fixed to 
Flat Earth Global Reference Frame 


6.3. Model Controllers and Issued Commands 
The commands in Table 6.1 can be issued to the 6DOF vehicle model. 


The command interface can be seen in Figure 6.8. 
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Table 6.1. 6DOF Model Commands 
































Command Range Units 
Propeller Speed 0 to 2000 RPM 
Depth -50 to 0 Meters 
Rudder Angle -13.6 to 13.6 | Degrees 
Forward HCTT Speed | -5000 to 5000 | RPM 
Aft HCTT Speed -5000 to 5000 | RPM 








REMUS Propeller RPM:Value 





PELUP PCCP Cee Pee eee e Pepe ee epee eee eee Peep eee pee ep eee tp eee ep eect peepee epee reper peenageens 
| 7 | | | | | | | | ! | | i nm at 


0 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 


Depth Command (Meters):Value 


aD A RC 


50 -48 -46 -44 -42 -40 -38 -36 -34 -32 -30 -28 -26 -24 -22 -20 -18 -16 -14 -12 -10 8 6 4 2 0 


Rudder Angle (Degrees):Value 





ERE CD CD CD DD 


13.6-12.6-11.6-10.6 -9.6 -8.6 -7.6 -6.6 -5.6 -4.6 -3.6 -2.6 -16 06 04 14 24 34 44 54 64 74 84 9.4 104114124 13.1 


FWD HCTT RPM:Value 


en Te 


5000 -4500 -4000 -3500 -3000 -2500 ‘2000 -1500 -1000 -500 0 500 1000 1500 2000 2500 3000 ‘3500 4000 4500 500 


AFT HCTT RPM:Value 





EC brad prrerprreey rout preerpererprresy 


5000 -4500 -4000 -3500 -3000 -2500 -2000 -1500 -1000 -500 0 500 1000 4500 2000 2500 3000 3500 4000 4500 500 


Figure 6.8. Command Interface for the 6DOF Hydrodynamic Model With 
Real Time User Specified Inputs (Set Here at 0 for All Parameters). 


These commands were selected because direct control of these parameters aids in the 
verification process. The only command in the 6DOF model that utilizes controllers that 


simulate physical controllers on the REMUS vehicle is the depth command. Since the depth 
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controller was verified with high confidence in the 3DOF model, the same controller was 
utilized in the 6DOF model. 


In an actual REMUS vehicle, the rudder angle and HCTT commands would be governed 
by a navigation/waypoint controller that would reduce the heading error of the vehicle 
attempting to achieve its next waypoint. This waypoint controller will need to be created 
and verified in further research. The RPM command in the physical REMUS vehicle is 
controlled by a speed controller. The scope of this research only includes the verification 
of the hydrodynamic model, but not the creation and verification of simulated controllers. 
To perform model verification, commands from a real world mission are sent to the DOF 
model, and the simulated vehicle’s behavior is compared to the actual vehicle’s behavior. 
The following experiments were used to verify the model’s behavior in response to these 


issued commands: 


1. Setting the rudder angle on the actual REMUS vehicle at various positions at a 
constant propeller RPM. The turn radius of the actual vehicle at various rudder angles 
and speeds is compared to the model’s turn radius. 

2. Importing real world REMUS propeller RPM and rudder angle data into the 6DOF 
model and comparing the vehicle’s and model’s speed. 

3. Setting various amounts of differential thrust on the forward and aft HCTTs on the 
actual REMUS vehicle while in a hover. The rotation/yaw rate is then recorded for 
the actual vehicle. This differential HCTT thrust will then be applied to the model 
and the rotation/yaw rate will be compared to the real world data. 

4. Real world vehicle propeller RPM, rudder angle, and depth commands are issued 
to the 6DOF model. The model’s pitch is then compared to the REMUS vehicle’s 


recorded pitch during the real wold mission. 


The results for these verification experiments are presented in Chapter 7. 


6.4 Complete 6DOF Model 

The complete 6DOF hydrodynamic model presented in Figures 6.9, 6.10, 6.11, 6.12, and 
6.13 represents a 6DOF variable speed hydrodynamic model of an AUV utilizing CTTs. 
This accurate model can be used to improve understanding of vehicle control for a wide 


variety of AUV operations in constrained environments, including terminal approach to a 
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docking station. Accurate 3DOF and 6DOF hydrodynamic models utilizing CTTs are an 
important capability for developing greater autonomy and enabling prolonged missions at 
sea. 
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Real World REMUS Commands 


Figure 6.9. Complete 6DOF Variable Speed REMUS Model Utilizing Cross- 
Tunnel Thrusters with Animation Subsystem, Real World Command Import 
Subsystem, REMUS Model's Real Time Control Subsystem. 


Figure 6.9 contains the REMUS Control Station (Figure 6.8), REMUS 6DOF Model(Figure 
6.7), and the REMUS Animator. The REMUS Animator takes the outputs from the 
REMUS 6DOF Model and creates a 3D representation of a REMUS vehicle, with a fly 
along camera angle, so that vehicle behavior can be critiqued and analyzed. The Real 
World REMUS Command block imports commands recorded during a real world REMUS 
mission and sends those commands to the model at the same rate the actual vehicle received 


the commands. 
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Figure 6.10. 6DOF Variable Speed REMUS Model Subsystems 


Figure 6.10 shows the subsystems contained in the REMUS 6DOF Model block, previously 
seen in Figure 6.9. The left block in Figure 6.10 takes the imported Real World Commands, 
or the user defined REMUS Control Station commands, and translates them into inputs that 
the 6DOF equations of motion can utilize. The 6DOF Model Force and Moment Equations 
Subsystem, represented in Figure 6.12, contains the equations of motion represented in 
Figures 6.1, 6.2, 6.3, 6.4, 6.5, and 6.6. The 6DOF Model Coordinate Transformation with 
Model Outputs Subsystem, shown in Figure 6.13, contains the model’s components that 
convert the force in body reference plane to position and motion in the flat earth reference 


frame. This transformation was shown in Figure 6.7. 
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Figure 6.11. 6DOF Model Controllers, Direct Force Equations, and Direct 
Moment Equations Subsystem. 
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Figure 6.12. 6DOF Model Force and Moment Equations Subsystem. 


71 


6DOF (Quaternion) 





Figure 6.13. 6DOF Model Coordinate Transformation with Model Outputs 
Subsystem 
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CHAPTER 7: 
6DOF Model Verification 





The 6DOF hydrodynamic model verification was conducted through four separate experi- 
ments. These experiments where designed to provide a reference of the REMUS vehicle’s 
behavior for comparison with the 6DOF model’s behavior. The purpose of these experi- 
ments was to provide the same command inputs to the 6DOF model as was seen by the 
real world REMUS. By comparing the model’s behavior with the REMUS vehicle’s actual 
response, the model’s fidelity was analyzed. These experiments are in no way a compre- 
hensive verification of the 6DOF model, however the results do provide insight into the 


capabilities and limitations of the created model. 


7.1 Surge Verification 

The surge verification experiment was conducted by providing the propeller and rudder 
angle commands from two REMUS missions to the 6DOF model. The vehicle’s speed, as 
recorded by the ADCP, was then compared to the model’s speed. The results are seen in 
Figures 7.1 and 7.3. 
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Figure 7.1. Speed of the 6DOF model vs. speed of the REMUS vehicle 
operating with CTTs and at various speeds 
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Figure 7.2. CTT RPMs During Surge Testing 
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Figure 7.3. Speed of the 6DOF model vs. speed of the REMUS vehicle 
operating without CTTs and at various speeds 


Figures 7.1 and 7.3 show that the mean speed of the model is approximately 4% to 9% higher 
or lower than the mean speed of the actual REMUS vehicle operating in Monterey Bay. One 
mission was conducted without CTTs activated and the other was conducted with CTTs 
activated. When considering that the currents seen by the REMUS vehicle were up to 4% 
of the vehicle’s overall speed, the behavior of the model’s speed is considerably accurate. A 
potential explanation for the speed discrepancy during the mission with CTTs enabled is that 
the CTTs created parasitic drag by increasing the turbulence around the REMUS vehicle, 
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possibly altering the coefficient of drag. The CTT RPMs during the CTT enabled mission 
is plotted in Figure 7.2. Both VCTTs and HCTTs were in operation at nearly maximum 
RPM, due to bang/bang control, for the entire mission. Constant CTT operation, along with 
altered fluid dynamic properties due to variations in the coefficient of drag, could explain 
why the REMUS vehicle’s mean speed is slightly less (approximately 9 %) than the (DOF 
model’s speed for a CTT enabled mission. More testing and experimentation is needed to 


confirm this hypothesis. 


7.2 Horizontal Cross-Tunnel Thruster/Yaw (Heading) Rate 


Verification 
To verify the behavior of the CTTs in the horizontal plane the REMUS vehicle was ordered 
to maintain a certain depth with little to no movement. This maneuver is called a hover. 
While the vehicle maintained its hover at the specified depth, a differential thrust command 
was sent to the horizontal CTTs. This differential thrust induced a rotation in the horizontal 
plane. These same depth and differential thrust commands were also given to the 6DOF 
model. The heading of the actual vehicle and the model were compared and the results are 


illustrated in Figure 7.4. 
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Figure 7.4. Heading of the 6DOF Model vs Heading of the REMUS Vehicle 
while performing a hover maneuver with differential commands being sent 
to the HCTTs 


The model’s yaw rate behavior is heavily influenced by the hydrodynamic coefficient N,,\;1. 


When N,,,| is properly estimated for the REMUS vehicle, an accurate yaw rate behavior is 
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witnessed, depicted in Figure 7.4. 


7.3 Pitch Verification 

The pitch of the vehicle during a mission is actively recorded. To verify the pitch behavior 
of the 6DOF model, the model was provided depth, propeller, and rudder commands from 
two REMUS missions in Monterey Bay. One mission was conducted without CTTs enabled 
and the other was conducted with CTTs enabled. Without CTTs, the REMUS vehicle’s 
pitch has only one actuator input, its dive fin, for controlling pitch. With CTTs enabled, 
both REMUS and the 6DOF model have more control authority over the vehicle’s pitch 
angle. When vertical CTTs are not activated, the vehicle has less control over its pitch. 
The pitch of the REMUS and the 6DOF model, over a 1000 second mission, can be seen in 
Figures 7.5 and 7.6. 
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Figure 7.5. Pitch of the 6DOF Model vs. Pitch of the REMUS Vehicle 
Operating with CTTs and at Various Speeds 
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Figure 7.6. Pitch of the 6DOF Model vs. Pitch of the REMUS Vehicle 
Operating Without CTTs and at Various Speeds 


The 6DOF model always had greater control authority over pitch than both real world 
REMUS missions. It is likely that this increased control authority over pitch is due to 
the assumption of benign environmental factors within the model and the assumptions 
made about vehicle configuration and trim. Nearly any environmental force acting in the 
Z-direction on the vehicle will induce a pitching moment increasing the pitch oscillation 
magnitude and frequency, especially if the force is variable. The model currently does 
not have the ability to model environmental factors, but environmental factors would have 


impacted the REMUS vehicle’s pitch behavior while it was operating in Monterey Bay. 


The assumptions about the modeled vehicle’s weight distribution/payload, and the behavior 
of the actual REMUS vehicle’s pitch controller (which is not modeled) are most likely the 
main causes of discrepancies in pitch behavior. The tuning of the real world REMUS’s 
pitch controller will heavily influence the magnitude and frequency of pitch oscillations, as 
well as the pitch overshoot. Pitch overshoot is clearly seen in Figure 7.5. The overshoot 
magnitude in Figure 7.5 is almost three times the value of the normal magnitude of pitch 
oscillation. A possible explanation for this large overshoot is poor tuning of the pitch 
controller in the REMUS vehicle with CTTs disabled. The vehicles owned by the CAVR 
most likely had their pitch controllers tuned for CTT enable operations. The poor pitch 
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behavior witnessed in the CTT disabled mission may be a result of the REMUS vehicle’s 


pitch controller not being tuned for a CTT disabled mission. 
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CHAP TER:¢: 
REMUS Controller Design 





A 3DOF non-linear dynamical model was presented in Chapter 3 and a 6DOF non-linear 
dynamical model was presented in Chapter 5. With the development of these new models, 
the design of controllers for the simulated REMUS vehicle is now possible. The controllers 
presented in this thesis were designed to be robust and adaptable, while also simulating the 


behavior of the physical controllers in the actual REMUS vehicle. 


8.1 Design Philosophy 

The controllers for the 3DOF and 6DOF models where split into two separate depth control 
autopilots, the dive fin deflection controller and the VCTTs’ controller (seen in Figures 8.1 
and 8.2). The pair of controllers are combined to provide the global vehicle depth control 
system, illustrated in Figure 8.3. The actual REMUS vehicle utilizes relatively simple field- 
tuned Proportional Integral Derivative (PID) controllers [7]. PID controllers were therefore 
utilized in the simulated vehicle. PID tuning was performed on the model until the behavior 
of the modeled REMUS closely approximated the characteristic dive and thrust behaviors 


of the actual vehicle. 
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Figure 8.1. Dive Fin Deflection Controller in REMUS Models 
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Figure 8.2. VCTTs’ Controller in REMUS Models 
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Figure 8.3. Combined Depth Controller in REMUS Models 


8.2 PID Control 

PID control is a classic control technique that provides a simple yet elegant method of 
controlling a dynamic system’s behavior through feedback. A PID controller reacts to the 
error between the input/command signal and a desired output signal. PID controllers have 
three elements, the proportional element (P), the integral element (I), and the derivative 


element (D). A PID controller can be represented by the following equation: 


Kp +; f edt + Ko (8.1) 


Where e is the tracking error for the system and Kp, K7, and Kp are the proportional, 


integral, and derivative gains of the controller. According to [7], Kp handles immediate 
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error and its magnitude impacts rise time and steady state error. K; impacts how well 
the controller can “learn” from past behavior. Increasing its magnitude can significantly 
reduce steady state error and lower rise time but can increase overshoot and introduce 
oscillations into the transient response of the system. Kp affects the system’s "anticipation" 
of future errors. Increasing the magnitude of Kp typically increases the system’s stability 
and dampens overshoot and oscillation during its transient response. Kp, K;, and Kp are all 
highly interactive, and altering one gain can have a significant impact on the efficiency of 
the other gains. While automated procedures for tuning PID controllers do exist and were 
attempted while tuning the controllers in the 3DOF and 6DOF models, the best results for 
modeled vehicle behavior were found when models’ controllers were field tuned to represent 


behaviors seen in actual REMUS missions. 


8.3. Tuning of the Simulated PID Controllers 

The model’s PID depth controller was tuned to produce simulated vehicle behavior that was 
comparable to the actual REMUS depth control performance. The tuning procedure was 
similar to field-tuning the PID gains on the actual REMUS vehicle. First, the proportional 
gain, Kp, was tuned to provide a similar response time as seen in the real world. Next, 
the integral gain, K7, was tuned to insure the proper offset between command depth and 
actual depth. Finally, the derivative gain, Kp, was tuned to remove any oscillations that the 


controller induced in the system’s behavior. 


Table 8.1. Depth Controller Gains Used in the 3DOF and 6DOF Models 
for Dive Fin Controller 

















Gain | Value 
Kp —800 
Ky 1 
Kp |2 
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Table 8.2. Depth Controller Gains Used in the 3DOF and 6DOF Models 
for VCTT Controller 

















Gain | Value 
Kp 800 
Ky 25 
Kp 15.1 














The values for the gains listed in Tables 8.1 and 8.2 were tuned and selected by constant 
experimentation and comparison between the models’ behavior and the REMUS vehicle’s 
behavior. Gains for the actual REMUS vehicle’s depth controller are field tuned in much 
the same manner in an effort to provide the desired vehicle behavior during its programmed 


missions. The REMUS vehicles’ controllers are field tuned by Kongsberg, the manufacturer. 


82 





CHAPTER 9: 
Behavior and Significance of the Coefficient of Drag 


(Cd) 





During the creation of the 3DOF model presented in this thesis, the Coefficient of Drag 
(Cd) significantly impacted the values of many other hydrodynamic coefficients. As the 
speed of a vehicle changes, Cd depends greatly on the size and speed of the vehicle in 
question. To create an accurate model of the REMUS vehicle, or any underwater vehicle, 
the value of Cd at various operating speeds must be known. Numerous calculations and 
verification through field testing in Monterey Bay, however, revealed that the Cd of the 
REMUS vehicle could be approximated as a constant value over almost the entire range of 
REMUS operating speeds. This realization significantly reduced the computational burden 
of a variable speed hydrodynamic model. This chapter provides a justification for this 
simplification, and suggests when this simplification is accurate (and appropriate) for use 


in the development of a variable speed model for an underwater vehicle. 


9.1 Behavior of the Coefficient 

It is known that the Coefficient of Drag, Cd, changes drastically with speed if the flow 
around the body is laminar. However, if the flow around the body of interest is turbulent, 
the coefficient of drag is relatively constant for nearly all speeds. This behavior can be seen 


in Figure 9.1. 
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Figure 9.1. Drag Coefficient Cd of Rough Circular Cylinders in Steady Inci- 
dent Flow for Different Surface Roughness Values k/D (k=Average Height 
of Surface Roughness, D=Cylinder Diameter, Rn = UoD/v, Us=Incident 
Flow Velocity.) Triangle, k/D = 110 * 107°; Circle, k/D = 450 * 10°; 
Square, k/D = 900 « 107°; Line. Source: [17]. 


Previous research on the REMUS vehicle found the Cd at a nominal operating velocity 
of approximately 1.5 m/s. This Cd was used in many single speed models and proved to 
be accurate at simulating the hydrodynamic behavior of the REMUS at around 1.5 m/s. 
However, to build an accurate variable speed hydrodynamic model of the REMUS vehicle 
that can simulate low speed docking, the Cd over an entire spectrum of operating speeds 
must be known. The first step in understanding the behavior of the Cd for the REMUS 
vehicle was determining the speed at which flow over the vehicle transitions from laminar 


to turbulent. Reynolds number can be calculated using equation 9.1, as described in [13]. 


p|U|D 
L 


Re= (9.1) 
After computing this parameter, Table 9.1 and Figure 9.1 can be used to determine the speed 
at which the flow on the REMUS vehicle is turbulent "enough" to assume Cd is constant. 
A Reynolds number of approximately 40,000 is the threshold where the REMUS vehicle is 


operating in the turbulent region. 
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Table 9.1. Behavior of Water in Certain Reynolds Number Regimes. Adapted 
from: [2], [13]. 





Re Range Behavior 

0<Re<l Highly viscous laminar creeping motion 

1 < Re < 100 Laminar, strong Reynolds number dependence 
100 < Re < 10° | Laminar, boundary layer theory useful 

10° < Re < 10* | Transition to turbulence 

107 < Re < 10° | Turbulent, moderate Reynolds number dependence 
10° < Re < co | Turbulent, slight Reynolds number dependence 



































The fluid interaction area has a significant impact on the Reynolds number. For the Reynolds 
number calculations seen in Chapter 2, the fluid interaction area was conservatively chosen 
to be the smallest area exposed in any body plane, which is the REMUS nose surface area 
that has a characteristic diameter of 0.190492 m. Even when using the most conservative 
fluid interaction area for the Reynolds number calculations, it was found that the flow over 
the REMUS vehicle remained turbulent until the vehicle’s speed dropped below 0.3 m/s. 


9.2 Suggestions for Model Development 
To find the critical velocity, |7|(m/s), of an AUV with a certain diameter D(meters) the 


following equation can be used, 


IT |(m/s) = [0.054757(m?/s)]/[D(meters)] (9.2) 
This equation was developed during this research by manipulating the coefficient of drag 


equation and assuming a Re of greater than 4000. The equation is provided for a quick 


reference for future researchers. 


The critical velocity for the REMUS vehicle with D=0.190492 m is, 


|T| = 0.288196m/s (9.3) 
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It is reasonable to assume that even while operating with CTTs, the REMUS is almost 
exclusively operating in the turbulent region. The experiments conducted with the 3DOF 
and 6DOF models, which utilized a constant Cd, and the data collected from the REMUS 
vehicle operating in the field, validated this assumption. Most torpedo shaped AUV’s 
will not be operating with CTTs. To maintain controllability of these AUVs, a minimum 
speed is required to ensure adequate fluid flow over their control surfaces. If a majority of 
torpedo shaped AUVs require at least a 0.5 m/s forward velocity to maintain controllability, 
the AUV diameter at which laminar flow is possible while operating at 0.5 m/s is about 
11 centimeters. Therefore, there is a strong argument to be made that the Cd for nearly 
all torpedo-shaped AUVs can be assumed to be constant across their range of operational 
speeds. This realization is valuable for any person creating or designing a hydrodynamic 
model of an AUV operating at various speeds. Also, more research is required to fully 
understand the impact of operating at very low speeds, where flow over the vehicle may be 


laminar and Cd can no longer be assumed as constant. 
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CHAPTER 10: 
Conclusion and Future Work 





10.1 Conclusion 


All models are wrong, but some are useful. 


-George E. P. Box 


This research was sparked by the problem of autonomously docking an AUV. To simulate and 
predict the behavior of a REMUS vehicle when reducing speed and utilizing CTTs, 3DOF 
and 6DOF hydrodynamic models were developed. These models were then programmed 
and analyzed using Simulink and were tested and verified against real world REMUS 


missions conducted in Monterey Bay. 


Both the 3DOF and 6DOF models provided "useful" models of the vehicle’s behavior. The 
assumptions made in these models do not impact their fidelity dramatically and allow for a 
simulated prediction of the REMUS vehicle’s behavior in the field. With the development of 
these models, a better understanding of the capabilities of a REMUS vehicle with CTTs has 
been accomplished. The models can simulate a REMUS vehicle’s maneuvering behavior 
while using CTTs during docking missions, simulate the vehicle’s depth control, simulate the 
vehicle’s maneuverability at various propeller RPM, and can be consulted when designing 
future REMUS docking stations. 


There were other unintended findings when developing the 3DOF and 6DOF models. The 
Cd was calculated and analyzed over the entire range of REMUS velocities. This study came 
to the conclusion that Cd of a REMUS vehicle can be assumed to be constant over nearly all 
the operational speeds. The creation of the models also provided the insight and catalyst for 
creating the pitch behavior parameter, Mpg. This parameter insures that the righting force 
acting in the pitch direction is always large enough to stabilize the vehicle’s pitch behavior. 
Observation in the field of the REMUS vehicle strongly supports the argument that there is 
a strong righting force due to the moment arm created by the center of gravity and center 


of buoyancy. 
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Unexplored hydrodynamic phenomenons were also uncovered during the verification pro- 
cess. The suspected reduction of vehicle velocity when operating CTTs due to parasitic 
drag is worthy of further study. More analysis and comparison between the model and the 
REMUS vehicle’s behavior can further justify that this phenomenon exists. Also tow tank 
analysis and computational fluid dynamics may provide a more concrete explanation for 
this parasitic drag witnessed in the initial comparisons between the model and the REMUS 
vehicle. 


10.2. Future Work 


The following goals can be pursued in future work: 


1. Calculation and verification of the minor coefficients used in the models need to be 
further analyzed. A more in-depth study of these coefficients may uncover further 
phenomenon and create a higher fidelity model of an AUV. 

2. Further development of a model that can accept environmental inputs such as wave 
action and currents will provide a needed capability for critiquing the feasibility 
of certain missions. Currents for a real world Monterey Bay mission are seen in 
Figure 10.1. The currents in Monterey Bay during the experiments were minimal, 
however their impact on the vehicle was still significant. Modeling and implementing 
environmental conditions would be a significant improvement in providing a higher 
fidelity and "useful" model. 

3. Revising the design and implementation of controllers to minimize the REMUS 
vehicle’s behavioral errors, such as the errors witnessed during the initialization 
period of the 3DOF model, would increase the model’s fidelity. A model that does 
not have this initialization period would provide a more accurate and faster prediction 
of the vehicle’s behavior, which would make it a better candidate for model based 
adaptive control. 

4. Integration of the DOF model into the REMUS vehicle’s secondary controller could 
provide a faster than real-time model to the controller, which may open more control 
options and provide an opportunity to test optimal trajectories prior to the execution 
of a maneuver. 

5. Development of a waypoint controller for the 6DOF model has already begun. Figure 
10.3 shows a mock waypoint mission, which utilized the DOF model developed in 
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this research interacting with a Latitude, Longitude, and Altitude controller, illustrated 
in Figure 10.2. This waypoint mission program was created to provide a top level 
mission simulation that accepted latitude, logitude, and depth commands and output 
overall vehicle behavior in the LLA reference frame. Model output is portrayed in 
Figures 10.4 and 10.5. This output allows for mission planning in LLA coordinates. 
The program provided a simulation of the vehicle’s trajectory in the mission area 
and an approximate time of completion for the mission. The actual REMUS also is 
programmed to achieve waypoints geolocated by latitude and longitude. Futhering 
this LLA capability and waypoint controller will provide useful tools for overall 


mission planning and asset allocation considerations. 


. Implementation of a neural network into the 6DOF model to learn and predict the 


REMUS vehicle’s behavior and to adjust coefficients for different payloads may 
provide a higher fidelity model for more than one vehicles’ configuration. 
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Figure 10.1. ADCP Estimated Currents During a REMUS Mission 
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Figure 10.2. LLA Waypoint Program For Enhanced Use Of 6DOF Model 
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Figure 10.3. Four Navigational LLA Waypoints Given To Enhanced 6DOF 
Model 
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Figure 10.5. Plot of the Vehicle's and Waypoint’s Latitude and Longitude 
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APPENDIX: Important Values 





A.1 Vehicle Parameters 


Table A.1. Vehicle Parameters 









































Vehicle Parameter Value Units 
Weight (W) 52.5 kg 
Buoyancy (B) 55.8362 N 
Length (L) 2.26 m 
Hull Frontal Area (A ,) 2.85 * 107? | m? 
Hull Projected Area (Ap) 0.430512 m 
Diameter of REMUS (D) 0.190492 m 
Fin Taper Ratio (¢) 0.654 

Fin Platform Area (S;in) 6.65 « 1073 | m? 
Moment Arm (X finpost) 3 —1.0475 m 
Fin Lift Slope (cia) sA2 

Pitch Behavior Parameter (Mpz) | 50 

















A.2. Environmental Variables 


Table A.2. Environmental Variables 





Environmental Parameter | Value Units 





Sea Water Density(p) 1030 kg/m 














Sea Water Viscosity (u)* | 0.00141 | kg/(m-s) 














‘Origin at Center of Buoyancy (CB) 
2At 10°C 
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A.3. Vehicle Hydrodynamic Coefficients 


Table A.3. 3DOF REMUS Model-Vehicle’s Hydrodynamic Coefficients 









































Vehicle Hydrodynamic Coefficients | Value Units 
Ky ~12.4759 | #4 
Xx -0.930 | kg 
Dx =98.6 (| 
Ziohil ~185.621 | *4 
Lies: -21.37 | 2, 
My —4.16 kg-m 
Muw 24.0 kg 
kg-m 
Mug -10.00 | 
Mow) 4.00357 | kg 
Mus, —22.3855 | 4% 
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Table A.4. 6DOF REMUS Model-Vehicle’s XYZ Hydrodynamic Coefficients 

























































































Vehicle Hydrodynamic Coefficients | Value Units 
ae, —12.4759 | 
Xi; -0.930 | kg 
xan -11.8 fg 
Xqq 71.93 a 
Xi 35.5 fg 
xe -1.93 tom 
V -28.6 | 8 
Yoo -2850 | 
Yer] 0.632 on 
Y, ~77.8 kg 
Y; 4.16 kgm 
ie 5.22 Aa, 
Vitis 9.64 aa 
Zi ~77.8 kg 
Zi 4.16 kgm 
Te =12.22 | #2 
Zalal -0.632 | 
Tei -28.6 | 
Zi ~185.621 | “2 
rae -21.37 | 2, 
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Table A.5. 6DOF REMUS Model-Vehicle’s MN Hydrodynamic Coefficients 










































































Vehicle Hydrodynamic Coefficients | Value Units 

kg-m 
My —30 rad 
Majal -188 | ap 
My —4.16 kg-m 
Muw 24.0 kg 

kg-m 
Mug —10.00 a 
My \w 3.18 kg 
Muué, —22.3855 | 42 
Nj 4.16 kg-m 
N; 4.88 ae 
Nw —24 kg 

kg-m 
Nur 2 gn 
Nolo} —3.18 kg 

2 

Nrtr| =245 | 
Nuué, —29,3855 | 42 
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